Josephson transport through a Hubbard impurity center. 
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We investigate the Josephson transport through a thin semiconductor barrier containing impurity 
centers with the on-site Hubbard interaction u of an arbitrary sign and strength. We find that in the 
case of the repulsive interaction the Josephson current changes sign with the temperature increase 
if the energy of the impurity level e (measured from the Fermi energy of superconductors) falls in 
the interval (— u,0). We predict strong temporal fluctuations of the current if only a few centers 
present within the junction. In the case of the attractive impurity potential (u < 0) and at low 
temperatures, the model is reduced to the effective two level Hamiltonian allowing thus a simple 
description of the nonstationary Josephson effect in terms of pair tunneling processes. 
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Josephson transport through a semiconductor barrier 
containing resonance impurity centers is a subject of in- 
tense current theoretical and experimental interest. Be- 
ing a coherent current of Cooper pairs, the Josephson 
current flowing through the impurity center is extremely 
sensitive to the presence of on-site impurity interaction 
[1,2]. This offers a possibihty of using the Josephson ef- 
fect as a unique spectroscopic tool for measurements of 
impurities energy states and calls for a theoretical study 
of Josephson transport through an impurity level with 
arbitrary strength and sign of the on-site interaction. 

The on-site Coulomb repulsion makes the occupation 
of an impurity center by a Cooper pair unfavorable. 
Thus, one would expect that the Josephson current flow- 
ing through an impurity level is considerably suppressed 
(unlike the normal current via a resonance state). How- 
ever, recent analysis of the hopping magnetoresistance 
data for different semiconductors, revealed a presence of 
a comfortable transmission channel via the double oc- 
cupied level (related to the states of the upper Hub- 
bard band) [3] with the very low, as compared to the 
naively expected value, repulsion energy. Such a repul- 
sion reduction may result from the polaronic effect, which 
sometimes can even "overscreen" the Coulomb repulsion 
reverting it to the effective attraction at the site. The 
so-called D'^'^-centers in semiconductors formed by sub- 
stitutional dopants in GaAs and AlGaAs alloys (see e.g. 
[4,5]) represent an example of those attractive impurities. 

In this Letter we investigate the Josephson transport 
through a semiconductor barrier with impurity centers 
of arbitrary on-site interaction strength and sign. We 
restrict ourselves to a sufflciently low impurity concen- 
tration such that the presence of an electron on a given 
site does not affect its neighbors. We consider the semi- 
conductor film to be thin enough, so that it's width is 
less than the average distance between the impurity cen- 
ters: in this case the transport is indeed determined by 
the tunneling through a single impurity rather than by 
hopping over a chain of resonant impurities levels. 

Given an on-site interaction strength, u, the occupancy 



of the impurity level is controlled by the level energy e. 
For u > (repulsion), the impurity level is double occu- 
pied if e < —u. It is single-occupied when —u < e < 0, 
and, finally, the impurity state is empty if £ > 0. In case 
of very strong repulsion u — > oo and zero temperature, 
the Josephson current changes sign abruptly as soon as 
the transition from the empty state to a single-occupied 
state occurs, in the later case a so-called 7r-junction is 
realized [1]. We show that this effect holds and that the 
jump of the Josephson current becomes even more pro- 
nounced for any finite positive interaction u. At finite 
temperatures in the regime —u < e < the double - and 
zero occupied states that carry positive Josephson cur- 
rents become exited. We demonstrate that it results in 
nonmonotonic dependence of the Josephson current on 
temperature: the current first increases from negative to 
positive values and then decreases. 

In case of an attracting center, the impurity level is ei- 
ther double occupied, if e < u/2, or empty when s > u/2. 
Close to the resonance, when e ~ u/2, the main contribu- 
tion to the Josephson current comes from pair tunneling 
processes, and the system can be described by the effec- 
tive two level Hamiltonian that includes pair tunneling 
processes only [2]. Solution of this model gives two en- 
ergy branches with energies E±{(f)) that depend on the 
phase difference (j) = (j)2 — 4>i of the superconductors and 
correspond to exactly opposite Josephson currents (see 
Fig. 3). At finite temperatures the upper level is exited 
and the Josephson current depends on the temperature as 
tanh£;(0)/2r with E{(j)) = S+((^) - S„(0). This model 
also allows for a simple enough analysis of the nonsta- 
tionary Josephson effect in the regime where the applied 
voltage V ^ min{A,u)/e. We will consider the simplest 
set-up where the Josephson current is shunted by the re- 
sistor and show that in this regime the current-voltage 
dependence exhibits resonant peaks. 

Model. We describe our system by the Hamiltonian 
of impurity level coupled with superconductors via weak 
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tunneling matrix elements ti and t2 

/c=l,2 

where Hi and H2 arc the BCS Hamiltonians of the 
superconductors and Hq is the impurity Hamiltonian 
^0 = £(^a da+U djdi . 

Perturbation theory in tunneling matrix elements. The 
current flowing through the impurity can be found as 
the expectation value of the current operator 7i = 
ti e i [-(/il^ (0) rfa — (iJjV-'ia(O)] flowing between the super- 
conductor 1 and the impurity. In the second order per- 
turbation theory in t\ the current I = is 



I' 



I=-Uetl\m / dTiFl {t)F4-t), 
Jo 



(1) 



where -Fij/>(t) is the local anomalous Green function of 
the superconductor 1, Fi^{u)) = ttz^i Ai/-\/|Aip + w^, 
where vi is the density of states. The anomalous 
Green function of the impurity center Jd(ri,T2) = 
{di{Ti)di{T2)) is found by the second order perturbation 
theory in t2 and the current I becomes 

I =-4ietltl j dr dra dn (t) 

xi^(0,r,r3,r4)F2(r4-T3), (2) 
where the correlation function 

K{Ti,T2,Ta,TA) = {d-\{Ti)dl{T2)d\{Tz)d\{T4^))d, (3) 

is defined with respect to impurity Hamiltonian Hq. As- 
suming for simphcity that |Ai| = IA2I = A we arrive at 
the final answer for the current in the form 

I = Io{ct>) {io + h e-^' + 12 e-^'-) Z-^ (4) 

where £2 = 2e + u is the energy of double-ocupied state, 

h{<l))=ATi'^e{tit2fviV2^H4'), (5) 

the impurity partition function Z = 1 + 2e~^^ + e~^^^, 
and the terms io and 12 are 

io = -A\e) + ]-B{e), 
£2 A 

i2 = —A^{£ - £2) + ]-B{£ - £2), 
£2 A 



with functions A and B defined by 
e A 



yl(e) = 2TV , , , 

B(e) = -AA2(e)/2e + ry „. 



(6) 
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FIG. 1. The .loscphson current /' = I/Io{(f>) as 
a function of the energy of the single-occupied impu- 
rity state e for u/A = 0.5 and different temperatures: 
T/A = 0, 0.02, 0.04, 0.1, 0.2 for plots 1-5 respectively. 
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Repulsive Hubbard center. In the case of repulsive 
interaction at zero temperature Eq.(4) simplifies to 
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£ < —U 
-U < £ < 
£ > 0, 



(8) 



that describes double-occupied, single occupied and un- 
occupied contributions respectively. The dependence of 
the Josephson current on e given by Eq.(8) is discon- 
tinuous (see Fig.l) with the jumps at points £ = and 
£ = —u where the Josephson current changes sign. The 
maximal values of the Josephson current Imax/Io{4') = 
\/u + (2/7r — 1/2) /A grow as 1/u when m ^ such 
that for u < t^v the perturbation theory in hopping ele- 
ments at £ « 0, —u becomes unapplicable. In the region 
— u < £ < the Josephson current docs not cxibit a sin- 
gular behavior because ocupation of the impurity level by 
one electron prevents from the resonant pair current flow 
even in case of small u. Thus the height of the jumps 
of the Josephson current at points £ = and e = —u 
decreases with the increase of the interaction strength u. 

At flnite temperatures the dependence of the Joseph- 
son current on e becomes continuous (Fig.l) and, even- 
tually, at high enough T the Josephson current becomes 
positive for all e. Thus for e laying within the interval 
(— u, 0) the Josephson current changes it's sign as temper- 
ature grows (Fig. 2a). This can be understood as follows: 
At low temperatures the current is given by the single 
occupied contribution ii which is negative, at finite tem- 
perature the zero and double occupied states correspond- 
ing to positive currents arc exited. Since the latter ones 
make larger contributions, the Josephson current even- 
tually changes sign with the increase of temperature. 



The contribution ii in (4) is 
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FIG. 2. a) The Josephson current /' = I/Io{4>) as a func- 
tion of temperature T' = T/A for different positive inteac- 
tions u and e — — it/2. b) Dependence of the Josephson cur- 
rent on e' = £ -|- u/2 for u/A = —0.5 and different tempera- 
tures: T' = 0, 0.05, 0.1, 0.2 (top to bottom). 



The temperature dependence in question implies the 
statistical averaging over impurity centers. If only few 
centers are present, one expects significant temporal fluc- 
tuations of the critical current. These fluctuations can 
be observed, for example, by including the Josephson 
junction in a superconducting loop. The energy of such 
a loop in the abscence of external magnetic field can 
be written as a function of magnetic fiux $ = LI, 
with L being the inductance of the loop, as Ei{^) = 
—Ic cos(27r<I)/$o) + *I'^/2i, where I,, is the critical current 
defined by Eqs.(4,5) and $o is the magnetic flux quan- 
tum. Assuming IcL » $o we see that for positive Ic the 
minimal solution corresponds to $ = while for Ic < 
the solution corresponds to $ = $o/2. Transitions be- 
tween different states of the impurity center lead to fluc- 
tuations of the magnetic flux through the loop. Change 
in the occupation number requires a quasiparticle tun- 
neling from the leads to the center or vice versa. Thus, if 
T < A the characteristic fliictuation time is of the order 
of r~-^ ~ {t^ /h)yexj£i{—A./T). The characteristic times 
spent in the empty, single-occupied and double-occupied 
states are given as tZ, rZexp(/3e) and rZexp(/3e2), re- 
spectively. 

Attractive Hubbard center. At u < the dependence 
of the Josephson current on the level energy e shows a 
resonance at values e « \u\/2, where the energy of the 
double-occupied state £2 approaches zero (see Fig. 2b). 
Near the resonance Eq. (4) can be simplified leaving only 
the terms that demonstrate singular behaviors 1/^2 



7 = 



Io{^) A\{e2 + \u\)/2) - A2((e2 - \u\)/2) e-^^- 



£2 1 e-^s2 

For T <C A the maximal value of the current is 

Imax=h{4>)A^{\u\/2)/2T. 



(9) 



This expression diverges at very low temperatures T < 
t^u where the approach based on the direct perturbation 
theory in tunneling elements becomes unapplicable. 



Effective low tem,perature model. In case of the attrac- 
tive interaction the Josephson current obtained by the 
perturbation theory in t\ and t^ has a form of the Gibbs 
average of two terms corresponding to unoccupied and 
double occupied states. This holds as long as T <C w, 
and the system in this regime can be described by the 
effective Hamiltonian [2,6] 

H^ff = e^b^b + ii [S+e*-^^ + be'''^'] + h pe''^^ + be'''''^], 

where b = did-^ is the hard-core boson operator satisfying 
[5, S'l'] = 1 - 2b'^b, b^ b^ = 0. The eigenfunction of the 
Hamiltonian Heff can be written as tp = a\0 > +/3\1 >, 
where |0 > is the empty state and |1 > is the double 
occupied (one pair) state. The operator of the current 
flowing through the impurity is 



i = 2eiii <ip\be 



Solving the Schrodingcr equation Hgf f ip 
two eigenstates with energies 



-bU''^'\ij>. (10) 

Eip we flnd 



E±(0) = £2/2 ±^S^ + El (cos(</.) + 1) (11) 

where 5"^ = e\l^ -\- (?i - ^2)^ El = 2*1*2. The currents 
corresponding to these two states are 



eEl sin((/>) 



^6^+El{cos{cl>) + l) 



Dependences of the energies E± and currents I± on the 
phase (j) are shown on Fig 3. Comparing the current J_ 
given by Eq.(12) with the result of perturbation theory 
in ti,t2 we find the effective tunneling matrix elements 
^1,2 = T^t\ 2^^1,2 ^(1^^1/2)- At T <C M the summation in (6) 
can be reduced to the integration resulting in 



A{e) = 



1 — (2/7r) arcsin(e/A) 
^l-e2/A2 



-A2 



< e < A 
e> A. 



(13) 



At finite temperatures the Josephson current is given 
by the thermodynamic average of the two states /(</)) = 
/_((/)) ta.n\\E((j)) /2T E{(j)) = E+{(j)) - E_{(j)). As in 
case of repulsive interaction, for the junction containing 
only few centers we expect strong temporal fluctuations 
of the current since the critical currents for two branches 
have opposite signs. In particular, for the weakly coupled 
superconducting loop controlled by a single center within 
the Josephson junction the flux fluctuations between the 
values $ = and $ = $o/2 are expected according to 
the discussion given above. 

Nonstationary regime. Let us consider the simplest 
experimental set-up (see Fig 4.): Superconductor con- 
tacts are shunted by the resistor while the potential of 
the superconductor 1 and the current Iq are fixed. The 
potential of the superconductor 2 V2 is determined by 
current conservation equation 
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FIG. 3. Impurity energy E' = (E — e2/2)/E{] as a function 
the phase = — 02 (black lines) and the corresponding 
Josephson currents I'j. = I±/{eEo) (dashed lines). 



V2/R + l2{t) = Io, 



(14) 



where l2{t) is the current flowing between the impurity 
center and the superconductor 2 

72 (t) = -2eit2lm{a*l3e-''^^} (15) 

The phase ^2 is related to the voltage V2 by 

<^2 = -2eF2. (16) 

To close the system of equations (14-16) we write the 
Schrodinger equation for the amplitudes a, /3 



ia = (3ti+l3t2 e"'*" 
i/3 = /3e. +ati+at2e''^^ 



(17) 
(18) 



The Cooper pair energy on the impurity site Si has a 
contribution from the potential of the superconductor 
2 = £2 + 2fceV2 where k = li/{h + h) and h^h 
arc the distances between the impurity center and su- 
perconductors 1 and 2 respectively. Eqs. (14-18) can be 
easily solved numerically. Solution for k = 0, e2/-Eo = 
0.5,?i/£'o = 2.0,12/ Eo = 0.5 and different shunting re- 
sistors are shown on Fig. 4 along with approximate adia- 
batic solutions shown by the dashed lines. The adiabatic 
approximation is valid in the low voltage regime where 
the phase (j)2 changes slowly and the current is given by 

hit) = -eEl sin(02t)/^<52 + E^{cos{(j>2t) + 1), (19) 

corresponding to the lower energy branch E^((p2)- At 
higher voltages the effect of Landau-Zenner tunneling 
between two energy branches cannot be neglected. As 
we sec it from the numerical solution, transitions be- 
tween the branches appear as resonances on the current- 
voltage dependence. Positions of these resonances can 
be easily estimated in case ti » 12] taking ^2 = in 
Eq.(ll) we obtain that the energy levels are splitted by 

AE = 2^e\/A + 1\. Resonances arise when the fre- 
quency of the phase oscillations of the superconductor 2 
u = 2eV2/h is related with AE as mhu = AE. Shown 
in the the inset of Fig. 4 is the differential resistance as 




FIG. 4. Voltage-current dependence of the shunted Joseph- 
son contact for receptivities Re^ /% = 1.6, 1.2, 0.8 (top to bot- 
tom). Upper inset shows the setup while the lower one shows 
the differential resistance (dV/dl ) as a function of the applied 
Voltage corresponding to the lower graph {Re^ /% = 0.8). 



a function of voltage. One clearly sees resonances satis- 
fying the above conditions. 

Conclusions. In conclusion we considered the Joseph- 
son transport through an impurity center for the cases of 
arbitrary sign and arbitrary strength of the Hubbard in- 
teraction. For repulsive centers we show that the Joseph- 
son current changes sign with temperature when the en- 
ergy of single occupied states e lays within the interval 
(— u, 0). If the junction contains only few centers strong 
temporal fluctuations of the current are predicted. In 
the case of attractive centers we study the nonstation- 
ary Josephson effect with the help of the effective model 
that takes into account only pair tunneling processes. 
We consider the case of resistivelly shunted Josephson 
junction and show that current- voltage characteristic has 
resonances associated with transition between two states 
formed due to coupling of the impurity with one of the 
superconductors . 
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